Magnetic structure of an imbalanced Fermi gas in an optical lattice 
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We analyze the repulsive fermionic Hubbard model on square and cubic lattices with spin im- 
balance and in the presence of a parabolic confinement. We analyze the magnetic structure as a 
function of the repulsive interaction strength and polarization. In the first part of the paper we 
perform unrestricted Hartree-Fock calculations for the 2D case and find that above a critical inter- 
action strength Uc the system turns ferromagnetic at the edge of the trap, in agreement with the 
ferromagnetic Stoner instability of a homogeneous system away from half-filling. For [/ < f/^ we 
find a canted antiferromagnetic structure in the Mott region in the center and a partially polarized 
compressible edge. The antiferromagnetic order in the Mott plateau is perpendicular to the direction 
of the imbalance. In this regime the same qualitative behavior is expected for 2D and 3D systems. 
In the second part of the paper we give a general discussion of magnetic structures above Uc- We 
argue that spin conservation leads to nontrivial textures, both in the ferromagnetic polarization at 
the edge and for the Neel order in the Mott plateau. We discuss differences in magnetic structures 
for 2D and 3D cases. 

PACS numbers: 03.75.Ss,67.85.-d,71.10.Fd 



I. INTRODUCTION 

Cold atoms constitute a promising route to simulate 
model Hamiltonians of strongly correlated many-body 
physics with accurate control of system parameter^l^. 
After major experimental breakthroughs with ultra- 
cold bosonic atoms like the Bose-Einstein condensa- 
tion (BEC) of alkali gase^^ or the observation of the 
superfluid-Mott insulator transition in a bosonic Hub- 
bard modeP, the field of ultracold atoms is currently 
addressing problems of strongly correlated fermionic 
systemaiSniS. Arguably the most prominent goal is the 
understanding of the phase-diagram of the fermionic 
Hubbard model, which is believed to be of major im- 
portance for high-temperature superconductivitjff"i2l A 
two-component Fermi-gas in an optical lattice is well- 
described by the single-band Hubbard model, whenever 
the energy gap to higher bands is much larger tha n on- 
site interaction, temperature, and chemical potentiaf32E!. 
Only recently the fer mionic Mott transition has been 
realized experimentalljU^'i^l xhe major challenge for 
studying magnetism of the fermionic Hubbard mo del is 
to reach temperatures below the Neel temperatur d^^ l ^° l 
In addition to the preparation of the antiferromagnetic 
state, characterization tools have to be developed to allow 
a clear identification of the magnetic structure. Po ssible 
experimental techniques include Bragg spectroscopji^^^, 
local meas urem ents of the magnetization^^ , noise 
correlation^^SHZ]^ p^. ^j^g recently realized quantum gas 
microscopd^H 

The experimental control of spin imbalance in Fermi 
gases offered a unique way to study pairing phenom- 
ena beyon d the standard BCS picture for attractive 
interaction^^Hm^ Motivated by these results, we address 
in this work the effect of spin im balance on the repulsive 
fermionic Hubbard modeI22123] While we study strong 



optical lattices, where a single-band Hubbard model is 
realized, the magnetic structure of weak to intermediate 
lattice strength including multiple bands has also been 
discussecP^. We find rich physics arising from the inter- 
play between antiferromagnetic and Stoner ferromagnetic 
instabilities and spin imbalance. 

The magnetic order of the two-dimensional repulsive 
Hubbard model has been extensively studied in the past, 
see e.g. ReffTUl Cold atoms in optical lattices differ 
in several ways from typical condensed matter systems. 
First, there is a superposed external confinement poten- 
tial, which divides the system in an incompressible Mott 
state in the center of the trap and a compressible region 
at the edge. Second, the total spin is conserved, which 
means that we need to minimize the energy of the system 
given a global magnetization rather than a finite Zeeman 
field. One interesting problem concerns the spatial distri- 
bution of the imbalance between Mott plateau and edge, 
and it turns out that the solution strongly depends on 
the interaction strength. The constraint of spin conser- 
vation affects the ferromagnetic instability at the edge by 
enforcing nontrivial spin textures'^^ which also affects 
the Neel order in the Mott plateau in the center, as we 
will discuss in section HVl 

In this work we study the repulsive fermionic Hub- 
bard model including a parabolic confinement poten- 
tial. In the first part of this work we perform unre- 
stricted Hartree-Fock calculations for the 2D case. Rel- 
evant physics for this system can be identified based on 
the mean-field phase diagram for the repulsive 2D ho- 
mogeneous Hubbard modeP^. Up to a critical inter- 
action strength Uc it predicts antiferromagnetic order 
close to half -filling and paramagnetic order elsewhere. 
In the spirit of a local density approximation one might 
then expect that cold fermionic atoms in an optical lat- 
tice have antiferromagnetic correlations in spatial regions 
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with one atom per site and are paramagnetic elsewhere. 
In order to account for a finite imbalance, the system 
has to change its magnetic structure. Using an unre- 
stricted Hartree-Fock approach for the 2D system we 
find a canted antiferromagnet in the Mott plateau in the 
trap center and a partially polarized edge. We note that 
canted antiferromagnetic order close to half-filling has 
been reported previously in Ref. 1321 With spin polar- 
ization along the z-direction, the canted antiferromag- 
net accommodates the imbalance forming a constant z- 
component of the local magnetization, and simultane- 
ously it benefits from the superexchange interaction by 
building up an alternating magnetic order perpendicu- 
lar to the z-direction. Fixing the global imbalance and 
increasing the interaction strength results in more imbal- 
ance flowing to the edge. 

Above a critical interaction strength IJc the unre- 
stricted Hartree-Fock calculation predicts that the sys- 
tem turns ferromagnetic at the edge of the trap, in agree- 
ment with the ferromagnetic Stoner instability of a ho- 
mogeneous system away from half-filling. Furthermore, 
the orientation of the antiferromagnetic order in the Mott 
plateau is perpendicular to the direction of the ferromag- 
net in the edge. Spin conservation has again a strong 
impact on the magnetic structure of the system, since 
a uniformly polarized ferromagnetic edge together with 
an antiferromagnetic Mott plateau are generally not al- 
lowed. We will discuss spin textures in 2D and 3D lattices 
for U > Uc, which fulfill spin conservation and which 
show the two prominent features predicted by the mean- 
field calculation, namely a) magnetic instabilities towards 
ferromagnetism in the compressible edge and antiferro- 
magnetism in the Mott plateau and b) at the interface 
between Mott plateau and compressible edge, the orien- 
tation of the antiferromagnet and the ferromagnet are 
perpendicular to each other. 

We are aware that the chosen mean-field approach gen- 
erally overestimates symmetry breaking and therefore the 
critical on-site interaction strength, Uc, corresponding to 
the appearance of an intrinsic ferromagnetic edge, will 
presumably be higher than the one predicted here. How- 
ever, intrinsic ferromagnetism away from half-fil ling is ex- 
pected for sufficiently large interaction strengthpSEZl and 
in fact experimental indications for itinerant ferromag- 
netism in a Fermi gas of ultracold atoms have been re- 
ported recently in Reff351 Given the tunability of the 
ratio between onsite interaction and nearest neighbor 
hopping, U/t, the interaction strength required for the 
presented phase separation should be accessible in ex- 
periment {U/t — 150 have been reported in Ref ITS)) . 

The paper is organized as follows. In Sec. [H] we intro- 
duce the model and in Sec. |III| we calculate the magnetic 
structure for U < Uc within an unrestricted Hartree-Fock 
approach. The topology of the intrinsically ferromagnetic 



edge arising for U > Uc is addressed in Sec|IV|and in the 
Appendix. Finally in Sec. |V] we summarize our findings 
and comment on the experimental significance of our re- 
sults. 



II. MODEL 

We consider the fermionic single band Hubbard model 
on a 2D and 3D cubic lattice with an external parabolic 
confining potential. The Hamiltonian is 



'4a '^jc 



(1) 



where a € {t, i} labels the two fermionic components, 
which are the eigenstates of the z-component of a spin 
algebra. These two components can either be the hyper- 
fine state of the trapped fermions or even correspond to 
different atomic species. Ci^ denotes the annihilation op- 
erator for a particle with spin a at site i, = cl^'^icr and 
Hi = Uia- are the spin resolved and total occupation 
of site i. U is the on-site interaction and t is the nearest 
neighbor hopping. Finally denotes the distance of site i 
from the trap center measured in units of the lattice spac- 
ing a and a = mw^a^/2 characterizes the strength of the 
external confinement. The associated energy scale is the 
confinement strength at the edge of the atom cloud with 
one atom per site, denoted by Vt . In 2D Vt — Na/n, 
where N is the particle number. 



III. UNRESTRICTED HARTREE-FOCK 
APPROACH IN 2D 

We now apply a Hartree-Fock mean-field decoupling 
in the spin and the density channel. Since the trap 
breaks translational invariancc, the mean-field parame- 
ters will be site-dependent. Allowing for arbitrary spin 
and density at each site we obtain the following mean- 
field HamiltoniarFll 



H — Ho + Hint 

Ho = -t cl^Cj^ + aY^ 



(2) 



where Si = {J2a,i3 4a^a.,i3Cii3) /2 denotes the spin oper- 
ator at site i ((? is the vector of Pauli matrices) and 
Mi = (Si) is the local magnetization. Magnetization 
and density are determined self-consistently for fixed to- 
tal particle number N. In the following we assume zero 
temperature. The energy of the self-consistent solution 
is given by the sum over the lowest N single-particle en- 
ergies of the Hamiltonian ^ plus the constant energy 

Eo = UEml-{n^)V^)■ 

An important subclass of self-consistent solutions are 
the ones with coUinear magnetization where My{i) — 
Mx{i) = on all sites. In particular, the generic phases of 
the homogeneous Hubbard model^ have a coUinear mag- 
netization; either ferromagnetic Mz{i) = M, antiferro- 
magnetic Mz{i) = (— 1)*M, or paramagnetic Mz{i) = 0. 
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However, we will show that generally the combination 
of trapping potential and imbalance will lead to a non- 
collinear magnetization profile. 

We are interested in the ground state for a given im- 
balance, characterized by the polarization P = (iV-|- — 
N^)/{N^ + N^), which is an experimentally controllable 
parameteil^. The imbalance is conserved, since the two 
components correspond to different internal states of the 
atoms (typically different hyperfine states) and transi- 
tions between these states are energetically forbidden un- 
less they are driven by additional lasers. The single par- 
ticle eigenstates of the Hamiltonian in Eq. ([2| only have 
well-defined spin if the magnetization is coUinear. Gen- 
erally, an expectation (5*2) 7^ can be tuned by spin- 
dependent chemical potentials or equivalently by a ficti- 
tious magnetic field in z-direction = —BSz- 

The parabolic confinement will decrease the density 
away from the trap center. In a local density approx- 
imation, a cross section through the trap corresponds 
to a cut through the (n, U) phase diagram at constant 
interaction U. Polarization can most easily be accom- 
modated by ferromagnetism, but also antiferromagnetic 
and paramagnetic regions can account for finite imbal- 
ance. As discussed in the introduction, in a canted 
antiferromaget a spatially constant component aligned 
with the field is added to the alternating component per- 
pendicular to the imbalance. The paramagnetic region 
can be partially polarized in the spirit of Pauli param- 
agnetism, where the polarization is proportional to the 
applied field. In the following we show that canted an- 
tiferromagnetic order is realized at half-filling, and we 
study how the imbalance is distributed between Mott 
plateau and edge as a function of interaction and im- 
balance. Self-consistent solutions of the Hubbard model 
on the two-dimensional square latt ice (l2| ) have either a 
coUinear or coplanar magnetizatiorPil^^and we can set 
My = without loss of generality. However, we note 
that enforcing vanishing global in-plane magnetization 
can lead to non-trivial three- dimensional topologies for 
the intrinsic ferromagnetP^l^ll^ which will be discussed in 
SecHVl 



A. The homogeneous system at half-flUing 

Fig. ([!]) shows the mean-field energies of canted and 
coUinear solutions as a function of increasing imbalance 
for the homogeneous system at half-filling. A rough ex- 
planation of why the canted antiferromagnetic order is 
favored can be given within the mean-field Heisenberg 
model. Here the energy increases only quadratically with 
polarization for the canted order, but linearly with po- 
larization for coUinear magnetization. Since the solutions 
are the same at the extreme values P=0 and P=l, the 
ground-state is always a canted antiferromagnet. 
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FIG. 1: (Color online) Energy per particle as a function of po- 
larization for the homogeneous Hubbard model at half-filling 
for various interaction strengths. Solutions with coUinear 
magnetization (Mz) have higher energy than solutions with 
canted antiferromagnetic (CAF) order. Results obtained on 
a 20x20 lattice with periodic boundary conditions. For the 
canted antiferromagnet, the imbalance is fixed by a fictitious 
magnetic field in the z-direction, however, the Zeeman energy 
is not included in the plotted energies. 



B. Magnetization profile in the trap 

Fig. ([2]) shows an example of a typical magnetization 
profile of a self-consistent solution at an intermediate in- 
teraction strength U = 5t. For the chosen parameters, 
the interaction is strong enough to form a Mott plateau 
with (n^i)) = 1 in the center. Furthermore, the trap 
strength, a — 0.02t, and the particle number, N — 540, 
correspond to Vt — 3At which is smaller than the on-site 
interaction so that double occupancies are absent. 

Within the Mott plateau we find canted antiferromag- 
netic order, as expected from the analysis of the homo- 
geneous system. The cross sections of the spin resolved 
densities and the local magnetization in panels (b) and 
(c) of Fig. ([2]) show that the edge is partially polarized 
and does not have antiferromagnetic order, although the 
x-component of the magnetization extends into the edge. 

We now consider the distribution of a fixed imbalance 
for various on-site repulsions. Fig. ([s]) illustrates that 
increasing interaction moves the imbalance to the edge, 
(we define the Mott plateau through jn^ — 1| < 0.05). 
Above a critical interaction strength (of order Uc ~ lOi) 
the edge is fully polarized and the Mott plateau is a pure 
antiferromagnet. The maximum in the majority density 
at the border of the Mott plateau can be understood by 
recalling that in the homogeneous system for strong in- 
teractions, there is a 1st order phase transition between 
an antiferromagnet close to half-filling and a ferromag- 
net at finite doping®. By decreasing interactions below 
Uc, the canting in the Mott plateau increases and the 
polarization at the edge decreases. 

Next we describe the magnetic structure as a function 
of the global polarization, P, keeping the other parame- 
ters fixed. For U = 5t, the upper panel of Fig. Q shows 
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FIG. 2: (Color online) Magnetic structure for U — 5t, 
P = 0.5, iV = 540, and a = 0.02 {Vt = 3.4t). (a) Local 
magnetization at the lattice sites. Vertical (horizontal) com- 
ponent of arrows encodes Mz (M^). Dark (blue) corresponds 
to negative values of Al^ and light (green) to positive values 
of Mx- X and y denote spatial coordinates, (b) Spin-resolved 
densities along cross section aX x = 0.5. c) Spin-components 
along cross section at x=0.5. The lattice size is 40x40. 



that both the polarization in the center with canted an- 
tiferromagnetic order and in the partially polarized edge 
increases linearly with the global polarization. The po- 
larization at the edge is always larger than in the center 
until the Mott plateau disappears close to full polariza- 
tion. 

We now discuss the case of strong interaction, i.e. 
U > Uc- Here the edge is intrinsically ferromagnetic. 
As shown in Fig. ([s]), at U ~ I2t the edge is already 
fully ferromagnetic in absence of any fictitious magnetic 
field that is otherwise used to fix a certain global imbal- 
ance. Given the total number of atoms in the trap, N, 
and the number of atoms in the edge, Nqi, this defines 
a critical polarization Pc = Nqi/N, which is Pc ~ 0.5 in 
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FIG. 3: (Color online) (a) Spatial distribution of imbalance as 
a function of interaction strength for constant global imbal- 
ance P — 0.5. Increasing interaction increases polarization at 
the edge (Poi) and decreases polarization in the center (Pi). 
The fictitious magnetic field applied to fix the imbalance is 
shown with stars, (b)-(e) Cross-sections at a; = 0.5 of spin- 
resolved densities and Mx- Labeling is shown in (e). 



Fig. ([3]) and Fig. Our mean-field approach predicts 
for P < Pc and U > Uc sl spatially uniform ferromagnetic 
edge with a direction other than the z-direction. This 
implies a finite global in-plane magnetization. However, 
as we discuss in detail in the next section such a solu- 
tion which is forbidden by spin conservation, and the 
preferred ferromagnetic order in the edge will have non- 
trivial spin textures for P < Pc and U > Uc- For now we 
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FIG. 4; (Color online) Spatial dependence of imbalance as a 
function of global imbalance. Poi denotes polarization at the 
edge. Pl (A''i) denotes polarization (number of atoms) in the 
Mott plateau for (a) U = 5t and (b) U = lOi. The dashed 
lines are explained in the text. All other parameters are as in 
Fig. ([2f. 



restrict our discussion to P > Pc and U > Uc- Then the 
ferromagnetic order at the edge points in the z-direction 
and the antiferromagnet in the Mott plateau is canted as 
shown in the lower panel of Fig [3] 

We now increase the number of particles so that the 
center of the trap is more than half-filled. In agreement 
with the symmetry of the homogeneous Hubbard model 
around half-filling, we find that the edge between the 
Mott plateau and double occupied sites shows similar fea- 
tures as the outer edge discussed above. Fig. ^ shows 
the magnetization profile and the spin-resolved densities. 
Here Vt = 15.7 which is larger than the chosen on-site 
interaction. The Mott plateau is formed on a ring and 
has canted antiferromagnetic order. Moving away from 
the Mott ring, the antiferromagnetic order rapidly van- 
ishes and the edge is strongly polarized. In fact, for this 
rather large value of U we see a small maximum of the 
majority component at the outer edge and a minimum 
in the minority component at the inner edge. 



IV. NON-TRIVIAL SPIN TEXTURES FOR 

u>u. 

The Hartree-Fock calculation predicts that above a 
critical interaction strength Uc the edge of the atom cloud 
turns ferromagnetic, even in absence of any imbalance or 
fictitious magnetic field. In the previous section we de- 
fined a critical polarization, P^, corresponding to a fully 
polarized ferromagnetic edge along the z-direction and 
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(b) Densities 
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FIG. 5: (Color online) Same as Figure [5] but for U = lOt, 
N = 2472, a = 0.02, and P = 0.37. The lattice size is 64x64. 



an antiferromagnetic Mott plateau. In this section we 
discuss qualitatively the magnetic structure for U > Uc 
and P < Pc. 

A cold atom experiment is prepared from a paramag- 
netic state with no optical lattice. Controlling the imbal- 
ance between the two fermion species, the initial state is 
characterized by 



(M,) = PN/2; (A4) = 0= (Mj,), 



(3) 



where P is the polarization and N the number of atoms. 
Since there is no coupling between the effective spin de- 
gree of freedom and the rest of the experimental sys- 
tem, the same constraints apply in the presence of an 
optical lattice and with strong onsite interaction [/IMIIII 
This additional constraint is always fulfilled in our mean- 
field treatment except for U > Uc and P < Pc, where a 
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spatially uniform ferromagnetic edge is predicted with a 
direction other than the ^-direction. However, such a 
solution leads to a finite global in-plane magnetization, 
which is forbidden by the boundary condition. In order 
to fulfill Eq. ([3]), itinerant ferromagnetism in cold atom 
systems can have non-trivial topology as shown recently 
for balanc ed sy stems with filling factor less than unity 
everywher^2lM]^ 

In the following we discuss the magnetic structure for 
U > Uc and P < Pc, both for two-dimensional (2D) and 
three-dimensional (3D) systems. We look for magnetic 
structures that fulfill spin conservation ([3| and which 
show the two prominent features predicted by the mean- 
field calculation, namely a) magnetic instabilities towards 
ferromagnetism in the compressible edge and antiferro- 
magnetism in the Mott plateau and b) at the interface 
between Mott plateau and compressible edge, the orien- 
tation of the antiferromagnet and the ferromagnet should 
be perpendicular to each other. Our qualitative analysis 
is based on the Ginzburg-Landau-type free energy func- 
tional (see e.g. Ref . l35l and 136]) 



E = 



^1 



(fr-\WM\ 



lAfl 



\Mo\y 



(4) 



where p is the positive stiffness constant, Mq is the mag- 
nitude of the favored magnetization, and /3 > deter- 
mines the cost of amplitude fluctuations. The favored 
spin texture for strong interactions, U > Uc, is deter- 
mined by minimizing the total energy under the con- 
straint of Eq. ([3|. In our qualitative analysis, we neglect 
that at the edge the system parameters in Eq. Q de- 
pend on the radius. This allows us to write the total 
energy of a spin structure, as a sum of three contribu- 
tions: the energies of the spin structures at the edge, 
inside the Mott plateau and at the interface of both re- 
gions. We note that the energy scale related with the 
spin structure of the ferromagnetic edge is of the order t 
and thus much bigger than the small superexchange t'^ /U 
that determines the spin structure in the Mott plateau. 
Therefore we first minimize the free energy of the intrin- 
sically ferromagnetic edge. The remaining two energy 
terms describe the interface between ferromagnetic and 
antiferromagnetic order at the edge of the Mott plateau 
and the energy of the spin structure in the Mott plateau. 
Based on the different scaling with the system size we ar- 
gue that the interface term dominates for large systems. 
While the interface term scales with r^~^, where tm is 
the radius of the Mott plateau and D denotes the dimen- 
sion, the antiferromagnet scales like In tm in 2D and like 
TMlnrM in 3D, as we show below. We minimize the in- 
terface term by choosing the orientation of the ferromag- 
netic and the antiferromagnetic order to be perpendicular 
to each other at the interface between Mott plateau and 
compressible edge. In the following we discuss solutions, 
where the Mott plateau has no net imbalance. In fact, in 
the limit of large interactions [/ — )• 00, the superexchange 
t^/t/ vanishes, so that one could allow for a strong polar- 
ization of the edge, by polarizing the Mott plateau in the 



opposite way. As estimated in Appendix [B] such a solu- 
tion is however higher in energy for realistic interaction 
strengths. 
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FIG. 6: (Color Online) Illustration of possible magnetic struc- 
tures ai U > Uc for balanced 2D systems with vanishing 
global magnetization. Dark (blue) arrows indicate spin tex- 
ture of the ferromagnetic edge. Light (green) arrows illustrate 
the magnetization in the antiferromagnetic Mott plateau. 
Black dots indicate the regular 2D lattice. As discussed in the 
text the energetically favored solution is depicted in (a) and 
consists of a vortex structure at the ferromagnetic edge (here 
in the xy-plane) that is perpendicular to the antiferromag- 
netic ordering in the Mott plateau (here in the z-direction) . 
(b) and (c) illustrate a Skyrmion and domain wall structure 
in the ferromagnetic edge. 



FIG. 7: (Color Online) Illustration of magnetic structures in 
2D for an imbalanced system with U > Uc and < P < P^. 
Outer (blue) arrows indicate spin texture of the ferromagnetic 
edge which consists of a vortex in the a^y-plane tilted towards 
the z-axis. Inner (red) arrows illustrate the orientation of 
the staggered magnetization in the antiferromagnetic Mott 
plateau which is perpendicular to the ferromagnetic order. 



We argue that (i) in presence of the Mott plateau a 
vortex structure for the ferromagnetic edge should be 
energetically favored as depicted in Fig [6] (a) , and (ii) a 
finite imbalance should result in a vortex structure of the 
ferromagnetic order parameter in the xy— plane together 
with a small z-component, see Fig. ([t]). An important ex- 
perimental consequence is a strong z-component of the 
antiferromagnetic order in the center, which is aligned 
perpendicular to the ferromagnetic order in the edge. In 
Appendix [A] we derive the energy of the different topo- 
logical orders of the ferromagnetic edge; vortex, domain 
wall, or Skyrmion. These structures are illustrated in 
Fig. ([g]). It turns out that for realistic parameters, the 
vortex is lowest in energy. For finite imbalance the edge 
will then be described by a ferromagnetic vortex in the 
xy— plane and a constant z-component. The energeti- 
cally preferred direction of the antiferromagnetic order 
in the Mott plateau is perpendicular to that of the fer- 
romagnet at the edge. The antiferromagnet in the Mott 
plateau will therefore have a small in-plane magnetiza- 
tion forming a vortex, which grows with increasing imbal- 
ance, and a strong z-component as illustrated in Fig. ([7| . 



B. 3D lattice 

Similar arguments can be applied to a 3D system. 
Taking into account the boundary condition of vanish- 
ing global magnetization in balanced systems and by 
applying Eq. Q, one finds that the preferred struc- 
ture of the ferromagnetic edge in a balanced system 
is a hedgeho g^^'"^ . As shown in Appendix |Aj the 
energetically preferred antiferromagnetic order in the 
center should then be either a planar vortex struc- 
ture with AIaf — ±Afoe0 or a 3D spherical vortex 



(b) 




FIG. 8: (Color Online) Illustration of possible magnetic struc- 
tures in 3D for a balanced system with U > Uc- Outer (blue) 
arrows indicate the magnetization at the ferromagnetic edge 
and inner (red) arrows illustrate the staggered magnetization 
in the antiferromagnetic Mott plateau. While the ferromag- 
netic edge always has a hedgehog structure the Mott plateau 
has either a planar vortex structure (a) or a 3D "spherical" 
vortex structure (b). While both structures have the same 
energy for the balanced system, the planar vortex is favored 
by finite imbalance, see Fig|9] 



Maf — MqCq, where = (— sin </>, cos 0, 0) and ee — 
(cos 9 cos (/), cos sin — sin 9) , are spherical unit vec- 
tors. 

Both solutions are illustrated in Fig. ([s]). They guar- 
antee that at the edge of the Mott plateau, where the 
antiferromagnetic order of the center of the trap has an 
interface with the ferromagnet order at the edge, the ori- 
entations of the antiferromagnet and the ferromagnet are 
perpendicular to each other. A violation of this require- 
ment would cost an energy that scales with the area of 
the interface r\j. Deformations of the perfect Neel order 
in the center of the trap, either in amplitude or phase, 
are minimized and the corresponding energy scales as 
tm ln(rAf /a). For perfectly balanced systems, the vortex 
within the Mott plateau could lie in any plane. Imbalance 
will deform the hedgehog leading to a net z-component, 
see Fig. While this does not affect the energy of a 
vortex in the plane it increases the energy for the vor- 
tices in other planes or for the spherical vortex. There- 
fore we expect that for imbalanced systems in 3D with 
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FIG. 9: (Color Online) Illustration of magnetic structures in 
3D for an imbalanced system with U > Uc and < P < P^. 
Outer (blue) arrows indicate the magnetization at the ferro- 
magnetic edge. Inner (red) arrows illustrate the staggered 
magnetization in the antiferromagnetic Mott plateau. Note 
that finite imbalance only deforms the hedgehog structure of 
the ferromagnet edge, while the antiferromagnetic order in 
the Mott plateau is unchanged, see Fig. [8](a). 



U > Uc, the antiferromagnetic order in the Mott plateau 
will form a planar vortex structure in the xy— plane as 
in Fig. ([9]). In contrast to the 2D case where we expect 
a strong z-component of the antiferromagnetic order for 
U > Uc and P < Pc, we expect a vanishing z-component 
in 3D. 



V. DISCUSSION 

In this work we studied an interacting two-component 
Fermi gas on a 2D and 3D cubic lattice subject to a 
parabolic external confinement. We analyzed the mag- 
netic structure as a function of the repulsive interaction 
strength and spin imbalance. Applying an unrestricted 
Hartree-Fock calculation for a 2D system, we identified 
the critical interaction strength Uc where the edge turns 
ferromagnetic and analyzed the spatial distribution of a 
finite imbalance between the two Fermi components for 
U < Uc- We found that the system has canted antiferro- 
magnetic structure at half-filling with antiferromagnetic 
ordering in the plane perpendicular to the imbalance and 
is partially polarized elsewhere. Fixing the global imbal- 
ance and increasing the interaction strength results in 
more imbalance flowing to the edge. We expect the same 
qualitative behavior for 3D in that regime. 

In the second part of the work we gave a general discus- 
sion of the magnetic structure above Uc both for 2D and 
3D. We showed that spin conservation generally leads to 
nontrivial spin textures, both in the Mott plateau and 
at the edge. We predict that the edge has non-vanishing 
in-plane magnetization with a vortex structure in 2D and 
a hedgehog structure in 3D. We furthermore expect that 



for U > Uc and small imbalance the antiferromagnetic 
order in the Mott plateau has a finite z-component in 
2D, while in 3D a vanishing z-component of the antifer- 
romagnetic order in the Mott plateau is predicted. 

We expect our findings to have clear experimental sig- 
natures if temperatures below the Neel temperature can 
be reached. A phase-contrast imag^^ showing the den- 
sity of each component separately can test our predic- 
tions of a Mott plateau with ferromagnetic borders. De- 
tection of a canted antiferromagnet in the Mott plateau 
requires direct access to the order parameter. This can be 
achieved for instance through noise correlation^^!! or by 
measuring the local magnetizatiorP^^. Additionally one 
can use Bragg spectroscopjEHHl where the double unit 
cell of the antiferromagnet results in additional Bragg 
peaks. Furthermore the intensity of the additional Bragg 
peaks can then be used to measure the strength of the 
z-component of the antiferromagnet. 
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Appendix A: Ginzburg-Landau theory 

Following Ref. 1351 we apply a Ginzburg-Landau type 
description of the magnetism based on Eq. (|4| to analyze 
the magnetic structure for U > Uc-, where the edge is in- 
trinsically ferromagnetic. By enforcing a vanishing global 
in-plane magnetization, the ferromagnetic edge acquires 
non-trivial topology. For the energy estimate we consider 
three energy contributions. The most relevant contribu- 
tion is the free energy of the intrinsically ferromagnetic 
edge. Thereafter the contribution of the interface be- 
tween ferromagnetic and antiferromagnetic order at the 
edge of the Mott plateau has to be taken into account, 
which is minimized by choosing the orientation of the fer- 
romagnet and the antiferromagnetic to be perpendicular 
to each other. Finally the free energy of the antiferro- 
magnetic Mott plateau has to be minimized. 

We simplify our calculation by assuming constant pa- 
rameters p, P, and Mo in Eq. thus neglecting a ra- 
dial dependence of these parameters due to the trapping 
potentiaP^.We denote the radius of the atom cloud by 
Rc and the radius of the Mott plateau by tm- 
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FIG. 10: Radial cross section through the ferromagnetic edge. 
Ferromagnetic edge starts at the border of the Mott plateau 
at r = vm and ends at r = Rc- While for the vortex the mag- 
netization does not change in radial direction, both Skyrmion 
and Domain wall do change in radial direction within a ring 
defined by ro < r < ro + L. Note that the radial component 
of the magnetization changes by 2n around the circumference 
while the z-direction is fixed. 



1. 2D lattice 

In a 2D system we expect the magnetization at the 
edge to form a vortex-like structure. Furthermore, we 
claim that for a small imbalance the vortex will lie in the 
xy-plane with a uniform magnetization component point- 
ing in the z-direction. The energetically preferred direc- 
tion of the antiferromagnetic order in the Mott plateau 
is perpendicular to that of the ferromagnct at the edge. 
At the interface, the antiferromagnet in the Mott plateau 
will have an in-plane magnetization forming a vortex and 
a z-component. The lowest energy corresponds to the 
maximally allowed z-component of the antiferromagnetic 
order parameter thus minimizing the in-plane vortex. 

We now give quantitative arguments for the physics 
described above based on a comparison of the energies 
of a ferromagnetic edge with different topologies; either 
a vortex, a domain wall, or a Skyrmion as depicted in 
Figs |6] and [TOj First we discuss the balanced system. 
For the vortex the direction of magnetization is indepen- 
dent of radius but it rotates by 27r on each circumfer- 
ence. A particular realization of a vortex is My = MoeV. 
However, for the balanced system there is global rota- 
tion invariance and the plane of the vortex is arbitrary. 
Using Eq. Q, the energy cost of a vortex is given by 
Ey — TTpM^n{Rc/rM)- Even in the absence of a Mott- 
plateau, the lattice spacing, ao, gives a natural cutoff for 
the core energy leading to Ey < TrpAfg ln{Rc/ao). A vor- 
tex naturally fulfills the requirement of vanishing global 
magnetization in all three spatial directions. 



Another possibility is the formation of a domain wall. 
In the inner ring rM < r < tq there is a uniform po- 
larization (e.g. M = MqCz) and within a finite region, 
ro < r < ro + L, the sign of the magnetization is inverted, 
e.g. M = Mo{l — 2(r — ro)/L)e*z. In the outer ring, 
ro + L < r < Rc, the magnetization points in opposite 
direction, e.g. M = —MqCz- While the inner and outer 
rings have a perfect uniform ferromagnetic order, the do- 
main wall is energetically costly due to the suppression of 
the amplitude of the order parameter. The energy cost 
is given hy Ep = np M§{ro/L + l/2)(4 + 4LV(15^2)), 
with ^ = ^ p/^PM"^) denoting the coherence length, tq 
and L are not independent of each other but related by 
the condition of vanishing global magnetization. In ab- 
sence of any Mott-plateau, tm = 0, the smallest allowed 
value is tq/L 0.6 which increases with tm- Neglecting 
the term containing the coherence length we therefore 
obtain a lower bound for the energy of the domain wall: 
Ed > TTpM^A. 

Finally, we estimate the energy of a Skyrmion. The 
magnetization is uniform (e.g. M = MqCz) in the inner 
ring, Tm < r < tq, and then it rotates by an angle an 
around a local axis in a ring of width L, tq < r < tq + L, 
e.g. M = Mocos(^^a7r)e; -I- Mosin(^^a7r)e;. For 
a = 1, the magnetization in the outer ring is inverted, 
while for other angles it has a vortex structure, e.g. 
M = Afo cos(a7r)e2 + Afo sin(a7r)£v. The Skyrmion in- 
terpolates between the inner and outer rings by tilting 
the order parameter, keeping the amplitude of the mag- 
netization fixed in constrast to the domain wall where the 
amplitude is suppressed. In the region, tq < r < Tq + L, 
the magnetization of the Skyrmion changes in the radial 
direction and along the circumference. The radial depen- 
dence of the magnetization gives rise to an energy con- 
tribution given by Es = npMQ{rQ/L + 1/2) (ctt^). Again 
the variables tq, L, and a are not independent of each 
other but related by the condition of vanishing global 
magnetization. By minimizing this energy only, and ne- 
glecting the energy cost of the change of magnetization 
along the circumference, we get a lower bound for the 
Skyrmion energy: Es > npM§{rM/{Rc - tm) + l/2)7r2. 

According to these estimates the lower bound of the en- 
ergy for the domain wall is larger than the total energy of 
the vortex for > exp(— 4)i?c ~ -Rc/50, and the lower 
bound for the Skyrmion is larger than the vortex energy 
for Tm > exp(— 5)ii!c ~ i?c/150. In fact, the real minima 
for both Skyrmion and domain wall will be larger. Since 
the radi us of the whole atomic cloud is about 50 lat- 
tice site^i^E^, our conservative estimate shows that the 
vortex should be favored for practically any size of the 
Mott-plateau. We note that for the results shown in the 
main part of this paper tm /Rc ~ 1/2- For the balanced 
system there is global rotation invariance and the plane 
of the vortex structure is arbitrary. However, in presence 
of a finite imbalance the energetically preferred magneti- 
zation will be a vortex in the zy-plane with a uniform fer- 
romagnetic z-component. Assuming that the directions 
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of the ferromagnetic edge and the antiferromagnet in the 
Mott-plateau are perpendicular, we expect the direction 
of antiferromagnet to have a large z-component. 



2. 3D lattice 

Minimizing the free energy in Eq. Q, one finds that 
the preferred structure of the fe rroma gnetic edge in a 
balanced 3D system is a 

hedgeho#SJ36| 

now explain 

why we expect the antiferromagnetic order in the center 
to have a planar vortex structure for U > Uc and P < P^. 
First, at the edge of the Mott plateau, the preferred di- 
rection of the antiferromagnetic order is perpendicular 
to the orientation of the ferromagnet order in the edge. 
A violation of this requirement will cost an energy that 
scales with the area of the interface r^j. In order to fulfill 
the boundary condition at the edge of the Mott plateau, 
the antiferromagnetic order in the trap center can neither 
have perfect Neel order nor a hedgehog configuration, 
since the latter needs to be oriented in the radial direc- 
tion. One possibility is to build up the 3D magnetization 
from the preferred 2D solution for each plane z, which 
is given by Maf = Mo {z/rMCp - [1 - (z/rM)^]^/^ej}, 
where Cz — (0,0,1) and Cp — (cos(0), sin((/)), 0) are 
cylindrical unit vectors. However, this solution is not 
realized in 3D, since the change in the z-component 
of the magnetization between different planes costs a 
large energy that scales with the volume of the Mott 
plateau Eaf ex. r\j/a?'. In fact, the preferred mag- 
netic orders in the Mott plateau are either planar vor- 
tex structures like Maf — Mqe^ or 3D solutions like 
Maf — MoCg, where — (— sin </), cos (/), 0) and eg — 
(cos 9 cos (p, cos 9 sin (p, — sin 9) are spherical unit vectors. 
For the balanced system these solutions have the same 
energy given by Eaf ~ 47rMoprAf ln(rM/a). However, 
imbalance will deform the hedgehog at the edge of the 
trap leading to a net z-component. Such a deformation 
increases the energy of these solutions except for a vortex 
in a;y— plane. Therefore we expect the magnetization pro- 
file in 3D for U > Uc and small imbalance to be given by 
a (slightly deformed) hedgehog ferromagnet at the edge 
of the trap and an antiferromagnetic order with a vor- 
tex structure in the xy— plane in the center of the trap. 
In contrast with the 2D case, where we expect a strong 
z-component in the antiferromagnetic order for U > Uc 
and P < Pc, we expect a vanishing z-component of the 
antiferromagnetic order in the Mott plateau in 3D. 



Appendix B: Polarizing the Mott plateau 

In section |IV| we propose spin structures for U > Uc 
and P < Pc that minimize the total energy while fulfilling 



the spin conservation The constraints ^ prohibits 
the formation of a uniform ferromagnet at the edge of the 
trap if simultaneously the Mott plateau has an antifer- 
romagnetic structure with zero net imbalance. However, 
since the constraints ^ apply to the whole system, one 
could imagine a system consisting of a fully polarized fer- 
romagnetic edge and a Mott plateau strongly polarized 
in the opposite direction, such that the global imbalance 
is small or even zero. We now justify why such solutions 
are energetically more costly than the ones proposed in 
section IIVI 



We therefore discuss the magnetic structure of a bal- 
anced Fermi gas in a 3D trap. The system can be divided 
into a Mott plateau for radius r < rjv/ and an edge for 
radius rjv/ < r < Re- In section |IV| we claimed that 
the ferromagnetic structure at the edge forms a hedge- 
hog. Applying the Ginzburg-Landau type free energy in 
Eq. Q, the energy of a hedgehog can be estimated. We 
strongly simplify our calculation by assuming a constant 
density at the edge.'^ The magnetitude Mq of the ferro- 
magnetic magnetization is therefore constant along with 
the stiffness, p, in Eq.Q. We now estimate the energy of 
the ferromagnetic hedgehog as Ep = Stt^Mq (i?c — tm). 
The stiffness of a homogeneous Fermi gas is given by 
p — l/(12fc|,xo) = h^/{36mn), where kp is the Fermi 
wavevector, xo the magnetic susceptibility, m the mass of 
the fermions, and n is the density (see e.g. Ref. |36). For 
sufficiently small densities the mass of a particle hopping 
between nearest neighbors in a 3D cubic lattice is given 
by TO = Tt? /{ta?')^ where a is the lattice constant and t the 
hopping matrix element between nearest neighbor sites. 
The stiffness on a 3D cubic lattice is therefore given by 
p K, to? /iiQn) and the energy of the balanced hedgehog 
becomes E « {2Tr/9)t{Rc — rA/)/(a'*n) ~ tN^/^ , where 
Np is the number of atoms at the edge of the trap. This 
energy could be gained by uniformly polarizing the edge. 
However, due to the conservation of the total imbalance 
the Mott plateau would then also be polarized by P^j ~ 
Np/Nm, where Nm denotes the number of atoms in the 
Mott plateau. The corresponding cost in energy can be 
estimated as Eaf « Plf2NM'U'^/U = 8t'^Nl/{UNM), 
where 2Nm is the number of nearest neighbors in the 
Mott plateau and 4£^ /U is the superexchange. The en- 
ergy cost of polarizing the Mott plateau is smaller than 
the energy gain of forming a uniform ferromagnetic edge 
if U/t > SGN^!^ /{ttNm). This is not satisfied for real- 
istic particle numbers Np,Nm > 10'^ and Np/Nm — 1- 
We thus conclude that in 3D for U > Uc and P < Pc 
the Mott plateau is not significantly polarized. The mag- 
netic structures that minimize the total energy and fulfill 
Eq. ^ are therefore the ones presented in section |I VI We 
expect similar arguments to hold in 2D. 
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